ITEP-LAT/2004-04 
MPI-2004-9 



A novel probe of the vacuum of the lattice 

gluodynamics 

M. I. Polikarpov^, S. N. Syritsyn^ and V. I. Zakharov* 

' Institute of Theoretical and Experimental Physics, B. Cheremushkinskaya 25, Moscow, 

117259, Russia 

JD \ * Max-Planck Institut fur Physik, Fohringer Ring 6, 80805, Miinchen, Germany 



O 
O 

(N 

X> 
IX, 



Abstract 

> 

We introduce a notion of minimal number of negative links on the lattice for a 
given original configuration of SU(2) fields. Negative links correspond to a large 
■ potential, not necessarily large action. The idea is that the minimal number of 

negative links is a gauge invariant notion. To check this hypothesis we measure 
correlator of two negative links, averaged over all the directions, as function of the 
distance between the links. The inverse correlation length coincides within the error 
bars with the lightest glueball mass. 



1 A traditional way to study spectrum of excitations is to measure correlators 



r—, ■ of various sources. For example, the correlator 1 : 

^ ; D gl (r) = (0|(G^(r)) 2 , (G^(0)) 2 |0) , (1) 

where G^ u is the gluonic field strength tensor, is sensitive to the glueball mass m g i 
at large (Euclidean) distances r: 

lim D q i(r) = const + (const) exp(— m„ir) . (2) 

r— >oo 

On the other hand, one can also study gauge dependent correlators, such as 

D% b (r) = (0|A«(r),4(0)|0) , (3) 

where A° is the gauge field. Such correlators are not unitary, generally speaking, 
and are not controlled by the glueball mass. In particular, the gluon propagator 
(jHJ) could even grow at large r. More realistically, i.e., as indicated by the lattice 
measurements, the gluon propagator falls off at large r but exhibits some spurious 
mass scales, for review see, e.g., (2]. For us, it is important that these spurious mass 
scales 9X6, cLS Si rule, lower than the glueball mass. 



1 For simplicity of notations, we do not account for the anomalous dimension. 
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2 In this note we introduce a new type of correlators which, as we hypothesize, 
might be unitary although they are not explicitly gauge invariant (like (fTJl) and 
check our hypothesis through lattice simulations of SU(2) gluodynamics. 

To explain the basic idea behind our measurements it is useful first to remind to 
the reader how one can introduce a gauge invariant condensate of dimension two in 
gauge theories One starts with the vacuum expectation value < (A®) 2 > which is 
obviously gauge dependent. One can, however, minimize this vacuum expectation 
value on the gauge orbits and the results < (A^) 2 > m in is gauge invariant by 
construction. To ensure that the minimum exists, the Euclidean signature is used 2 . 

We generalize this idea to the case of a Z(2) projection of the original SU{2) 
fields. In this projection the link variables take the values ±1 (for review see, e.g., 
[3]). Upon the Z{2) projection there is still remaining Z(2) gauge invariance. We 
fix this gauge freedom by minimizing the number of negative links over the whole 
lattice. We speculate, furthermore, that the density of these negative links might 
well be gauge invariant, in analogy with < {A a ^) 2 > m in- 

Moreover, we expect that the correlator of the negative links is unitary, like (|T]). 
This speculation goes apparently beyond the ideas discussed so far. Indeed, the 
negative links in the continuum limit correspond to singular potentials, A a ~ 1/a, 
where A is the lattice spacing. Naively, one could argue that such fields are artifacts 
of the lattice and, as a manifestation of this, their correlator dies off on distances 
of order a. Why we call these arguments naive, is because the so called P-vortices 
are defined as unification of all the negative plaquettes in the Z(2) projection and 
exhibit sensitivity to the scale of Aqcd, for review see jSJ. It is in analogy with 
these observations that we expect the "propagator of negative links" to scale in the 
physical units as well. This is of course a daring possibility which cannot be proven 
a priory but only supported or rejected by measurements. 

4 In more detail, we perform measurements both in the Direct- and Indirect- 
Maximal Center Projections (DMCP and IMCP). The details of calculations are 
given in the Appendix. As a result of SU(2) —* Z(2) projection, the original SU (2) 
field configurations get projected into the closest configuration of Z{2) gauge fields. 
The remaining Z(2) gauge freedom is then fixed by maximizing the functional 

F{Z) = Y J Z x , ll (4) 

X,fl 

with respect to Z(2) gauge transformations {Z x ^ — > z x Z x ^ z x +^, z x = ±1). 

After gauge fixing only the positions of the negative links are relevant and it is 
reasonable to change the variables: 

Z x ,n = {1, if %x,fi = -1; if Z x>fl = 1} (5) 

Moreover, to imitate a scalar correlator on the discrete variables, we consider here 
the isotropic correlator defined as 

G^v(r) = — <Z 0fl Z x ^>, (6) 

T r<\x\<r+§ 

2 A modification suitable for the Minkowski space is to work in the Hamiltonian formalism 0]. 
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where the summation is over all links Z x ^ for x lying in the spherical layer r < 
\x\ < r + |; N r is the total number of links in this layer. 

The correlator (jSJ) tends to ci noil- vanishing constant £r(oo) — <C Z x ^ u > 2 as 
r — > oo (< > is the average density of the negative links) and we fit the data 
by the expression (see (J2J)): 

G(r) = G(oo) + Cexp{-mr}. (7) 

We thus get the mass parameter m for various values of the lattice spacing a. 
Logarithmic plots for the isotropic correlator G(r) are presented in Fig. \$\ior IMCP, 
while the corresponding values of the mass m are depicted in Fig. [21 for IMCP and 
DMCP. 

5 As is seen from Figs. 11121 the correlator G Atl/ (r) scales in the physical units 
(within the error bars) and the mass parameter is close to the scalar ++ glueball 
mass (m(0 ++ ) = 1.65 db 0.05 Get" [7]). Thus, our measurements support the idea 
that the correlator © is in fact gauge invariant and unitary. Of course the results 
are not analytical but pure numerical and, in principle, the picture can change at 
smaller values of the lattice spacing. The non-triviality of this observation is that 
it is a correlator of potentials, which are not explicitly gauge invariant. Moreover, 
in the continuum limit the negative links correspond to singular potentials. 
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Appendix 

We perform our calculations both in the Direct jS] - and the Indirect (Hj - Maximal 
Center Projections (DMCP and IMCP). The DMCP in SU(2) lattice gauge theory 
is defined by the maximization of the functional 

F 1 (C/) = ^(TrC/^) 2 , (8) 

with respect to gauge transformations, U x » is the lattice gauge field. Maximization 
of JHJ) fixes the gauge up to Z(2) gauge transformations and the corresponding Z(2) 
gauge field is defined as: Z x> u = signTrC/^. To get IMCP we first fix the maximally 
Abelian gauge maximizing the functional 

F 2 (U) = J2^{U x ,^ 3 U+^ 3 ) , (9) 

with respect to gauge transformations. We project gauge degrees of freedom U(l) —* 
Z(2) by the procedure completely analogous to the DMCP case, that is we maximize 
Fi(U) (JBJ) with respect to U(l) gauge transformations. 
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Table 1: Parameters of configurations. 



p 


Size 


NlMCP 


Ndmcp 


2.35 


16 4 


20 


20 


2.40 


24 4 


50 


20 


2.45 


24 4 


20 


20 


2.50 


24 4 


50 


20 


2.55 


28 4 


37 


17 


2.60 


28 4 


50 


20 



To fix the maximally Abelian and direct maximal center gauge we create 20 
randomly gauge transformed copies of the gauge field configuration and apply the 
Simulated Annealing algorithm to fix gauges. We use in calculations that copy which 
correspond to the maximal value of the gauge fixing functional. To fix the indirect 
maximal center gauge from configuration fixed to maximally Abelian gauge and to 
fix the Z(2) degrees of freedom one gauge copy is enough to work with our accuracy. 
We work at various lattice spacings to check the existence of the continuum limit. 
The parameters of our gauge field configurations are listed in Table D To fix the 
physical scale we use the string tension in lattice units jHj, \fo = 440 MeV. 
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Figure 1: ln( 
for IMCP. 
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1) vs lattice spacing ioi ji = v (a) and for /i 7^ 1/ (b), results obtained 
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Figure 2: Mass parameters for IMCP and DMCP for jx = v (a) and for jx ^ v. 
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